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ABSTRACT 

We find restrictions on the mass and mixing angle of the dark matter sterile neutrinos 
using X-ray observations of Coma and Virgo galaxy clusters. 



1 INTRODUCTION 



Recent experiments on neutrino oscillations (for a re- 
view see e.g. Strumia & Vissani (2005)) determine non- 
zero difference of masses of various species of neutrinos. 
The experiments on atmospheric neutrinos (Fuknda et al. 
1998) are explained via — > i/t oscillations with mass 
difference Arriatm = [2.2to;|] • 10"^ eV^ while experi- 
ments (Ahmad ct al. 2002; Eguchi ct al. 2003) determine 



mass squared difference Am^oi — [8.2 



,+0.3l 

-o.sJ 



lO""" eV^ for 



J^e — > I'/i.T transitions. Probably, a simplest way to consis- 
tently explain these data is to add to the Standard Model 
several (at least two) gauge singlet fermions - right handed, 
or sterile neutrinos. However, the absolute values of masses 
of active and sterile neutrinos are not fixed by these exper- 
iments. 

It was noticed long time ago (c.f. Dodclson & Widrow 
(1994)) that a sterile neutrino with the mass in the keV 
range appears to be a viable "warm" DM candidate. If such 
a neutrino is a main ingredient of the DM, it is potentially 
detectable in various X-ray observations. To wit, there exists 
a (sub-dominant) radiative decay channel of sterile neutrino 
A'^ into active neutrinos Va and photon with energy E = 
ms/2, with the width (Pal & Wolfcnstcin 1982; Bargcr ct al. 
1995) 
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Here is the mixing angle defined as 
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with A'Id being the Dirac mass and rris is the sterile neutrino 
mass. 

Although even for the main decay channel (A^ 
3ua) the lifetime of sterile neutrino must exceed the age 
of the universe, the decay A^ Va'y can result in a 



potentially detectable X-ray flux (Dolgov & Hansen 2002; 
Abazajian ct al. 2001), since the density of DM in the uni- 
verse is large. For example, the contribution to the diffuse 
X-ray background (XRB) can be comparable with measured 
XRB and is given by (see Boyarsky ct al. (2005) and refer- 
ences therein) 
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where Pqj,,//o are the dark matter density in the universe 
and the Hubble constant. Neutrinos decaying at different red 
shifts produce a broad X-ray line with extended "red" tail. 
Such a feature in the XRB spectrum is, in principle, readily 
detectable (and distinguishable) from the broad-band con- 
tinuum of observed XRB (Grubcr ct al. 1999). The non- 
detection of the DM decay feature in the XRB signal en- 
ables to put a bound on 9 and rUs roughly at the level 
of (Boyarsky et al. 2005) 



Q.sSLn{2e) < 3 X 10" 



1 keV 



(4) 



where Q,a is the present day density of sterile neutrino, un- 
derstood as DM candidates. 

Clustering of the DM at small red shifts results in the 
enhancement of the DM decay signal in the direction of large 
mass concentrations, such as galaxy clusters. Taking into ac- 
count the typical overdensity TZ = p/puM of a galaxy cluster 
is at the level of 7?. ~ 10^, while typical cluster size is about 
D ~ Mpc ~ W-''^Ho\ one can flnd that the DM decay flux 
from a galaxy cluster. 



8nDl 



(5) 



: npl^,DDlQ^°^ is the mass of DM within telescope's 



field of view (FoV) and Dl , Dg are the luminosity and an- 
gular diameter distances to the cluster) is comparable to the 
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background DM decay signal 



■ TIDHq ~ 1 



(6) 



Although the total DM decay flux from the XRB and from 
the cluster are of the same order, their spectra are different. 
The flux from the cluster would be detected as a narrow 
line whose width AE — AEdat is determined by the spec- 
tral resolution of an X-ray detector. At the same time, the 
DM decay contribution into XRB is produced by the de- 
cays at red shifts 2 ~ 1 and, as a result the DM decay 
line is broadened to AE ~ ma/2. Thus, in spite of the fact 
that the compact DM sources at 2; ~ give just moder- 
ate enhancement of the DM decay flux, the enhancement of 
the signal in the narrow energy band centered on the line 
energy E = ms/2 could be large. From the energy depen- 
dence of dFxRB /dE (see e.g. Boyarsky et al. (2005)) one 
can estimate that the enhancement can be by a factor of 
ma/ AEdct ~ 10 — 100 for instruments of spectral resolu- 
tion of Chandra or XMM (in this paper we will use XMM 
observations of Virgo and Coma clusters) . 

However, the detection of the DM decay line in the 
galaxy clusters is complicated by the fact that most of the 
galaxy clusters show strong continuum and line emission 
from the hot intracluster gas. The typical temperatures of 
the intraculser gas. 



G TZpD^iD^nip ~ 10 
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1 Mpc 
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lie exactly the range of Chandra or XMM-Newton satellites 
- 1-10 keV. In fact, for a typical nearby X-ray bright cluster, 
the diffuse continuum emission from the cluster core is by 
a factor ~ 10^ — 10'^ stronger than the XRB emission from 
behind the core. Thus, the factor 10—100 increase of the flux 
in the DM decay line from the cluster direction is "washed 
out" by the strong increase of the X-ray continuum emission 
from the direction of a cluster. We will see however than in 
spite of such a strong background, the bound on sin^ 29 can 
be improved by a factor of about 2-4 in comparison with 
XRB (4). 

Previously this question was addressed in the paper 
by Abazajian ct al. (2001); Abazajian (2005). There the 
authors claimed the following limit on the neutrino mass 
and mixing angle, derived from the non-observation of the 
DM decay line from the center of the Virgo galaxy cluster: 
sin^(26l) < 1.0 X 10"^ ( keV/ms)"* for 1 < < 10 keV. This 
limit is several times better than Eq. (4). 

In this paper we re-analyze the bounds on the sterile 
neutrino mass and mixing angle imposed by the observa- 
tions of Coma and Virgo galaxy clusters. As we will discuss 
in Section 4, our constraints are weaker than those claimed 
by Abazajian ot al. (2001); Abazajian (2005). We will com- 
pare in details these two results in Section 5 and explain the 
origin of the difference. 

In this work we show that the problem of the large 
continuum X-ray emission from the cluster core can be 
partially relaxed, if one searches for the DM decay signal 
from a region outside the bright X-ray core of the clus- 
ter. The reason for that is the following: outside the core 
of the cluster, the surface brightness proflle of intraclus- 
ter medium (ICM) behaves as 5(r) ~ r^" (where, e.g. 
for Virgo cluster a ~ 1.6 (Fabricant & Gorcnstcin 1983; 



Schindlcr ct al. 1999; Young et al. 2002)), while the sur- 
face brightness profile of the DM is more shallow 5dm (j") ~ 
(see e.g. C'avalicrc & Fnsco-Fcmiano (1976), for mass 
density of Virgo cluster). Therefore, by moving away from 
the center, one can put much tighter constraints on the 
DM parameters. For the case of Coma cluster observa- 
tions with XMM one needs to use periphery observations 
("coma mosaic" observations, performed in May- June 2000, 
(Bricl ot al. 2001; Neumann ot al. 2003)). For the case of 
Virgo this goal can be achieved already for the XMM obser- 
vation of M87 (observation ID 0114120101) Indeed, as the 
core radius of Virgo cluster is very small (ro < 1.6'), the 
existing observations with XMM EPN camera (field-of-view 
~ 15') allows one to study properties of Virgo cluster till 
about lOro. We also show that the bound can be further 
improved if the X-ray spectra of galaxy clusters are ana- 
lyzed by the method proposed by Boyarsky ot al. (2005) for 
the search of the DM decay line in the XRB spectrum. We 
describe the analysis of the data from Coma and Virgo ob- 
servations in Section 3 and present results in Section 4. 



2 DM DECAY LINE VS. THERMAL 
BREMSSTRAHLUNG EMISSION. 

The main obstacle for improving the bound on the neutrino 
parameters from the X-ray background by using galaxy clus- 
ters observations is strong X-ray continuum and line emis- 
sion from the cluster core. In this Section we show that this 
obstacle could be partially avoided if one extracts the X-ray 
spectrum from the regions beyond the core of the continuum 
X-ray emission. The idea is that the radial surface bright- 
ness profile of the continuum emission is normally steeper 
than the radial surface brightness profile expected from the 
DM decay line. 

Indeed, the surface brightness profile of a re- 
laxed cluster of galaxies is usually well-fitted by the 
so-called isothermal /3-model (Cavalicro & Fusco-Fcmiano 
1976; Sarazin & Bahcall 1977): 

So 



5(r) 



(8) 



(1 + (r/ro)2)3-3-o.5 ' 

where r is the projected distance from the center of the clus- 
ter, <So is the overall normalization factor and ro is the "core 
radius" of the region of X-ray emission. Parameters So, ro 
and P are different for different clusters (see Section 3 be- 
low). 

The volume emissivity of the intracluster gas which re- 
sults in the surface brightness proflle (8) is 

Vo 



V(r) 



(9) 



(1 + (r/ro)2)3/3 

where r is the (3d) distance to the center of the cluster. 
The continuum X-ray emission from the intracluster gas is 
the thermal bremsstrahlung for which the volume emissivity 
is proportional to the square of the gas density. The radial 
density profile of the hot intracluster gas is therefore 

no 



(10) 



(1 + (r/ro)2)3/3/2 

To calculate the radial mass profile M(r) (gas plus galax- 
ies plus the DM) of the cluster one has to make an as- 
sumption that the intracluster gas is in hydrostatic equilib- 
rium (Cavalicro & Fusco-Fcmiano 1976; Sarazin & Bahcall 
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1977). One can really expect this assumption to hold only for 
relaxed galaxy clusters, but as it turns out even for e.g. Virgo 
this model gives good predictions (c.f. Nulscn & Bohringer 
(lOf)."))). Under this assumption the overall mass profile can 
be calculated from the Newton's law 

- '''''''^'^'^^ + ^ ^"^^ dr - 7^ ' 

(11) 

where T(r) is the radial temperature profile of the hot intr- 
acluster gas. In the region where the temperature does not 
change significantly, one can calculate the overall mass pro- 
file analytically by substituting (10) into (11). The assump- 
tion of hydrostatic equilibrium them immediately leads to 
the following total mass dependence on the distance r from 
the center of the cluster (Cavalioro & Fusco-Fcmiano 1976; 
Sarazin t Bahcall 1977) 



M(r) = 1.13 X 10"Mq (i- 



{r/rof 



1 keV Mpc 1 + (r/ro)2 
The overall density profile is given by 

3 + {r/rof 



Ptotal(T') — PO, total 



(l + (r/ro)2y 



(12) 



(13) 



In the first approximation one can estimate the DM density 
as Pdm = Kptotai with K < 1. 

The surface brightness profile (8) should be compared 
with the surface brightness of the DM. Integrating the vol- 
ume emissivity of the DM decay line over the line of sight 
one can find the brightness profile of the DM decay line 



5dm (r) = 



(14) 



where M-DM[r) — 2 f^ dz UumWr'^ + z^) is the column den- 
sity of the DM as a function of the distance from the cluster 
center. Substituting ptotai from (13) one obtains 



5dm (r) 



2 + [r/raf 



8(1 + (r/ro)2)3/2 



(15) 



Integrating 27rr-5DM(T') (15) from r = to some r one finds 
the fiux of DM from a FoV circle with the projected radius r: 



iW = 4—^2 Iro) 



(16) 



where we have defined a geometry factor g{r), describing 
that part of the cluster, which we are using for the mea- 
surement. In case of the circular region with the projected 
radius r one gets from (15): 



(r/ro) 



(17) 



g{r/ro) = 

One can see that 5dm (r) ~ 
profile of DM decay line is significantly more shallow than 
the profile of the continuum X-ray emission from the hot 
intracluster gas: 



for r ^ ro. Thus, the radial 
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r » ro . 
(18) 

As it is clear from Eq. (8), /3 < |, therefore, the share of DM 
contribution to the total brightness (at energy E = ms/2) 
increases as one moves away from the center of a cluster. 



3 DATA ANALYSIS 

In this work we have chosen to analyze data of XMM- 
Newton observations of Virgo galaxy cluster (observation 
ID 0114120101, June 2000) and Coma cluster. For the 
latter we have taken two observations which constituted 
the part of the so called Coma mosaic (Briel et al. 2001; 
Neumann ct al. 2003): that of the center of Coma clus- 
ter (observation ID 0124711401, May 2000) and one of the 
peripheral observations ("Coma 3", ID 0124710301, June 
2000). 



3.1 "Total flux" restrictions. 

First, we find the restrictions on the parameters of sterile 
neutrino, based on the total emission from the ICM. Namely, 
using the XMM observation data and estimating a mass of 
DM in the center of a cluster as in Section 2, we restrict 
parameters rris and sm^(26) by demanding that the flux of 
DM decay line (16) did not exceed the total flux in the 
energy bin, equal to the 3cr where with of the Gaussian line 
a is equal to the spectral resolution of XMM. We will often 
dub this method simply total flux restrictions. This method 
provides the most robust exclusion. 

3.1.1 Analysis for Coma cluster 

To find the typical numerical values of the DM decay fiux as 
compared to the thermal bremsstrahlung fiux, let us use the 
parameter values for the Coma galajcy cluster. Its surface 
brightness profile is well defined by the /3-model (8) with 
/3 = 0.75 and ro = 10.5' = 0.3 Mpc^ (see e.g. Bricl ct al. 
(1992); Neumann ct al. (2003)). The virial mass of Coma 
cluster, contained in the virial radius rcoma = 3.5 Mpc is 
Mc'oma = 1-3 X lO^'^ M© (Briel et al. 1992). Substituting 
these numbers into Eqs. (12)-(13) and taking fraction of 
gas in the Coma cluster 18% (Briel ct al. 1992) one finds 
that 



po.DM ~ 3 X 10 



Mq 

Mpc3 



(19) 



Substituting these values in Eq. (16) one finds that the fiux 
of DM from a FoV circular region with projected radius r is 
given by 



erg 



(20) 



FDM.Coma ^ 6.7 X 10"* 0' [t^I 9ir)- 

L 1 keV J cm^ ■ s 

where g{r) is given by Eq.(17). Flux (20) should be com- 
pared with the emission of the ICM in the Coma center 
region. We use observation of Coma center with EPN cam- 
era of XMM- Newton satellite (observational ID 0124711401) 
to extract flux per energy bin, equal to the spectral resolu- 
tion of XMM (ABdet ~ 200 eV) in the interval 0.5-10 keV. 
The geometry factor for the circle of 700" is equal to 0.82. 
The resulting restriction on nis and sm^{2d) is shown on 
Fig. 1 in the red solid line. Qualitatively this result can 
be described as follows. The flux (per energy bin) from 
Coma center^ at 1 keV -Fi keV = 1-62 x 10"^^ erg/cm^-s 

^ We take h = 0.71 and the distance to the Coma cluster Dj^ = 
98 Mpc (red-shift z = 0.023). 

2 Compare Bricl ct al. (1992); Fi keV ^ 1-1 x lO'^erg/cm^ ■ s in 
200 eV energy bin. 
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Figure 1. Exclusion plot from tlie total flux restriction from the 
circular region with radius 700" in the center of Coma cluster. 
The straight green lines represent sin^{29)m^ = const lines, cor- 
responding to the total flux at energies 1 keV and 10 kcV (in the 
200 eV energy bin). 
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Figure 2. Exclusion plot based on the restrictions derived from 
a total flux from Coma 3 peripheral observation (red solid line) 
as compared with Coma center (green line, identical to the red 
curve on Fig. 1). Thinner blue lines represent sin^(26) = const 
with normalization corresponding to fluxes at 1 keV and 7 keV. 



is about an order of magnitude bigger, than that at 10 keV 
Fio keV = 2.5 X 10"^^ erg/cm^ • s. The flux of DM Eq. (20) 
is bounded by these two fluxes and exclusion curve in the 
space (tos, sin 29) is bounded by the two straight lines 
rnf sin^(20) = const (dashed lines on Fig. 1), where const is 
defined (via Eq.(20)) by the fluxes Fi kov and Fio kcV corre- 
spondingly. 



3.1.2 Coma periphery 

As discussed in Section 2, to improve the bound on total 
flux, we move away from a center of the cluster. Unfortu- 
nately, modern X-ray telescopes are normally "narrow FoV" 
instruments and, e.g. the XMM pointing toward the center 
of Coma cluster covers only the region of the size r ~ ro. 
Therefore we use one of the observations from Coma mosaic 
(May- June 2000) (Briol ct al. 2001; Neumann ot al. 2003). 
Although the flux from Coma 3 periphery falls by a fac- 
tor 20, as compared to the Coma center'^, we see that rela- 
tive change between 1 keV and 10 keV is again an order of 
magnitude (c.f. Fig. 3). For Coma 3 observation (distance 
between the centers of two observations is 39' and FoV ra- 
dius is ~ 12'), the geometry factor in Eq. (20) is equal to 
ffComas ~ 0.18 which improve an exclusion plot Fig. 2 by 
about factor of 5. 



3.1.3 Analysis for Virgo cluster 

We repeat the analysis of the previous Section for the 
case of center of Virgo cluster (M87). For this cluster 
again the surface brightness is well-fitted by the /3-model 
with /3 ~ 0.4 and core 
and So = 5 X IQ-^^ 

198.3; \'ouug ct al, 2002). The mass in the center of the 



cm^ -s-arcmiii 



radius ro — 1.6' = 8 kpc 
^ (Fabricant & Gorenstcin 



Virgo cluster can be described by the isothermal /3-model 



For Coma 3 observation Fi koV - 
Fio kcV ~ 10"^^ erg/cm2 ■ s. 
" We adopt the distance to M87 Dl 



2 X 10 ^'^ erg/cm^ 
18 Mpc {z = 0.0044) 
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Figure 3. Comparison between Coma center (red curve) and 
Coma 3 (green curve) fluxes (per 200 eV energy bin). 



(see e.g. Nulscn & Bohringer (1995)). Using various mea- 
surements of the mass distribution in the center of the 
Virgo cluster (Nulscn & Bohringer 1995; Evrard et al. 1996; 
McLaughlin 1999), we obtain 

,16 Mo 



po,DM ~ 6 X 10 



Mpc^ 



(21) 



and correspondingly flux from the FoV circle with (pro- 
jected) radius r is given by 



7.7 X 10" 



erg 



(22) 



The geometry factor g{r) is defined in (17). 

Using the M87 observation by XMM-Newton (ID 
0114120101), we extract the total fiux from the central 11' 
of the EPN FoV (centered at M87 galaxy).'^ The geometry 
factor for this case is gMg7(ll') = 6.8. 

Using expression (22) for the flux of DM we get a 

^ The actual experimental data is shown on Fig. 9 in red line 
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Figure 4. Restrictions based on the total flux from the circle 
with the radius 11' with the center at M87 (red line). Green line 
represents the data for Coma periphery - same as on Fig. 2. Blue 
line - restriction (4) from Boyarsky ct al. (2005). 



restriction on parameters of sterile neutrino. The thermal 
bremsstrahlung of ICM around M87 has temperature about 
T = 2.54 keV, which is about 4 times smaller than that 
for the Coma cluster. Therefore, the flux drops by about 
two orders of magnitude between 1 keV and 10 keV.® This 
leads to an exclusion plot, which is weaker than that from 
a Coma periphery at small energies, but gets stronger, than 
the former restriction at masses nis > 12 keV, as shown on 
Fig. 47 



3.1.4 Virgo periphery 

As we said above, in general periphery of the cluster does 
not fit into one XMM FoV. However, Virgo cluster is a spe- 
cial case in this respect. Indeed, in spite of the fact that 
Virgo cluster is the closest to us galaxy cluster, the an- 
gular size of its core is rather small by the cluster scales. 
The surface brightness of the core of the Virgo cluster ra- 
dius ro = 1.6' Fabricant & Gorcnstein (1983); Young et al. 
(2002) and thus the FoV of the XMM EPN camera covers 
about lOro. Nevertheless, the ratio of 5gas/5DM is changing 
rather slowly - as (r/ro)""'^* 

On the other hand, taking ring (with the radii 9' and 
11'), rather than the full circle with the radius ~ 11' leads to 
the decrease of the geometric factor in Eq. (22) by a factor 
of about 5, which compensates the gain almost entirely and 
therefore, moving away from the center does not improve 
the restriction, based on the total flux from the vicinity of 
M87. 



^ The measured flux at 1 keV is Fi i;ev — 4.8 X 10~^^ erg/cm-^ ■ s. 
This number should be compared with the one, obtained from 
the surface brightness profile of Fabricant <fe Gorcnstein (1983), 
F ~ 4.3 X 10-" erg/cm2 ■ s. 

One should take all total flux restrictions at the energies above 
bremsstrahlung cut-off with the grain of salt. Indeed, as will be 
discussed in more details in Section 3.2.1-3.2.2, the uncertainties 
of measurement of this flux is comparable with the flux itself. 



3.2 Statistical analysis of the XMM data 

In order to further improve the bounds on the sterile neu- 
trino parameters (as compared to Sections 3.1), we apply 
the method suggested in Boyarsky ct al. (2005). Namely, we 
use the fact that the spectral shape of the DM decay line 
may differ significantly from that of the ICM continuum and 
therefore adding a line which is "too strong" (e.g. with the 
flux equal to that of the continuum) will simply contradict 
observations. 

In general, typical X-ray spectrum of a galaxy clus- 
ter is dominated by the bremsstrahlung continuum and 
line emission from the hot intracluster gas. To obtain 
a good fit of observed data to physical model (such as 
e.g. MEKAL (Mcwc ct al. 1986; Liodahl ct al. 1995) or 
Raymond-Smith (Raymond & Smith 1977)), one needs to 
take into account that their parameters (e.g. temperature 
and abundances) vary with radius and even with azimuthal 
angle; that the ICM plasma is often described by multi- 
component models; etc. Therefore, often the sophisticated 
data analysis is required to extract the "best-fitting" pa- 
rameters for these models and one can only fit these models 
to the data, coming from small annular sectors of the FoV. 

However, for our purposes one needs any "best-fit" 
model to the data, gathered from previously described annu- 
lar or circular regions of FoV. We are not interested in val- 
ues or distributions of temperatures, abundances, etc. of the 
ICM plasma. Rather, we want to put restrictions on param- 
eters of the sterile neutrino, based on non-detection of the 
presence of the decay line against an emission background. 
Both from experimental data and from theoretical models, 
one can see that X-ray spectrum of galaxies has a form of 
continuum plus separate lines. Knowledge of physics (the 
correct choice of model) allows to relate various parameters 
of the model (such as positions and intensities of the lines, 
abundances, etc.) to each other, however, the structure of 
the model remains the same: continuum emission plus lines. 

Therefore, we introduce a phenomenological model for 
the cluster spectrum which provides a good fit to the X-ray 
data, but at the same time has a limited number of param- 
eters so that one can easily constrain the maximal possible 
contribution of the DM decay line into the X-ray emission 
from the cluster. The restriction is much stronger in those 
regions of the spectrum, where one has no lines - the non- 
detection of the decay line against monotonic continuum, 
essentially puts the restriction on the flux of the line com- 
parable to the error of the continuum flux. At the regions, 
where emissions lines are present, the restriction of statis- 
tical analysis is much weaker and the flux of DM decay is 
essentially equal to the excess flux above the continuum. 

The models and method of their construction are de- 
scribed in Appendices A-C. These models provide a satis- 
factory fit to the observed X-ray spectra of Coma and Virgo 
clusters. Obviously, addition of a strong line feature with the 
integrated line flux equal to the model flux within a given 
energy bin would destroy the fit to the data. Moreover, it 
would provide a non-satisfactory fit already at much smaller 
DM line flux, of the order of the statistical error of the data 
in the narrow energy bin. 

In order to quantify the above qualitative observation 
we use the following algorithm to flnd the ma^ximal possible 
contribution of the DM decay line at a given energy: 
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Figure 5. Red line: exclusion plot based on statistical analysis 
from Coma center observation (circular region of the radius 700" 
around the center of the cluster). Green line - restriction based 
on the total flux from the same region (curve from Fig. 1). 



Figure 6. Red line: exclusion plot based on statistical analysis 
from Coma 3 peripheral observation of Section 3.1.2. Green line 
- restriction based on the total flux from the same region (curve 
from Fig. 2). 



(i) We fix the parameters of a phenomenological model 
to their best fit values. 

(ii) We introduce an additional model component which 
is a Gaussian line at the energy Eq with the width AE equal 
to the energy resolution of the instrument (70 — 150 eV for 
XMM EPN camera and 15%(£/60 keV)"" '' for Beppo-Sax 
PDS). 

(iii) Then we increase the normalization of the additional 
Gaussian until it starts to distort the overall fit to the data. 
As a simple criterium for this we use the increase of reduced 
X^, which leads to the vanishingly small null hypothesis prob- 
ability. The actual value of Ax^ depends of course on the 
number of degrees of freedom and has to be chosen individ- 
ually for every data set. We confirm by the visual inspection 
that for such a distortion, fitted model does not follows the 
data anymore in the region of E ^ Eq. 

(iv) We make a "scan" over the line energy 0.5 keV ^ 
Eo ^ 10 keV (and 15 keV < 40 keV for Beppo- 
Sax data) to find the maximal possible normalization of the 
Gaussian as a function of energy. 

As discussed above and in Boyarsky ct al. (2005), the only 
assumption about the empirical model, which we make, 
is that the power-law regions of the spectrum are actual 
monotonic bremsstrahlung and do not contain "dips" see 
Fig. Alb. 

In the Sections 3.2.1- 3.2.2 below we compare this 
method with that of the Section 3.1 for the case of Coma 
and Virgo cluster. 

3.2.1 Bounds from statistical analysis of the Coma cluster 
observations. 

The simple model for Coma center X-ray emission is pro- 
vided in Appendix A. The results of statistical analysis (and 
demonstration of its improvement over the total fiux restric- 
tions) are shown on the Fig. 5. One can see that overall 
one gets about an order of magnitude improvement over 
the results of Section 3.1.1. One can also see "spikes" at 
E ~ 1 keV, E ~ 6.5 keV, etc. {m, ~ 2 keV,13 keV,etc.) - 
positions of the strong emission lines in the Coma spectrum 



(see Appendix A for their identification). This is expected. 
Adding a line in the energy region where data is monotonic 
and well described by the smooth power-law bremsstrahlung 
emission distorts the fit much more than adding the line with 
the same strength to the region of multiple element emission 
lines. 

Similarly, the results of statistical analysis of the data 
from Coma 3 region are shown on Fig. 6. Here, the improve- 
ment is not so drastic as in case of Fig. 5. The reason for that 
is obvious - the statistics for Coma 3 is lower and therefore 
the fitted curve can be distorted stronger, without spoiling 
reduced x too much. From looking at the Fig. 6 one can see 
that the proposed statistical analysis is in some sense com- 
plimentary to that of the total flux analysis: namely, while it 
is stronger at small masses (rUs < 10 keV, it gets more and 
more weak as increases - tendency which is quite oppo- 
site to that of the total fiux restriction. The reason for that 
is clear - at large energies there is almost no X-ray emission 
from cluster, and the difference between a very weak signal 
from X-ray emission and background has very large uncer- 
tainty. Therefore, total flux restriction for rUs > 10 keV is 
provided only for illustrative purposes and the results of the 
statistical analysis should be considered as a correct limit.* 

3.2.2 Bounds from statistical analysis of the of Virgo 
cluster. 

Let us apply the same strategy to the XMM pointing to- 
ward M87/ Virgo cluster. For this we consider the spectrum 
collected from a thin ring 9' < r < 11' around the M87. We 
attempt to add a DM decay line to the resulting spectrum, 
described in Appendix C and find the maximal strength of 
the line that does not deform the spectral fit (in the sense 
described in the previous sections). The resulting exclusion 
plot is shown by the dashed lines in Fig. 7. We see that the 

* We have used the Beppo-Sax observation of Coma cluster, to 
obtain constraints in the range 15 keV ^ E ^ AO keV. The ob- 
tained restrictions (both statistical and "total flux" ) are about an 
order of magnitude weaker than restrictions from XRB. 
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Figure 7. Restrictions from tlio statistical analysis of the pe- 
riphery (ring 9'— 11') around the M87 (green dashed line). This 
restriction is about an order of magnitude better that the one, 
coming from the total flux, measured from the same region (red 
solid line) at small energies {mg < 8 keV and then become worse 
(red dashed line), as statistical error of the data grows. This is 
due to the fact that for i? > 4 keV total flux is measured with 
uncertainty comparable or exceeding the flux itself. The straight 
line represent the XRB restrictions (4). 
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Figure 8. Exclusion plot for parameters of sterile neutrino, based 
on analysis of the cluster of the Coma and Virgo clusters of galax- 
ies. The red solid line represents the restriction, coming from the 
Coma periphery ("Coma 3") observation, while the green dashed 
line, coming from observation of Virgo cluster, provides stronger 
restriction for rrie > 14keV. The blue dashed curve represents the 
result of statistical analysis of the Coma center. The straight ma- 
genta line represents the XRB restrictions from Boyarsky ct al. 
(2005). 



statistical analysis improves the restrictions over that of the 
total flux for small masses (i.e. for which E = ms/2 < T, 
where T ~ 2.5 keV is the temperature of the ICM for the 
center of the Virgo cluster, while gets significantly worse for 
higher energies. As discussed in Section 3.2.1, this is due to 
the fact that the errors in flux are comparable (or even ex- 
ceed) the flux itself and therefore the results of the statistical 
analysis should be considered as the correct limit. 

One can also see that the restriction from the Virgo is 
not as strong as the one, which we obtained from the Coma 
center (shown on Fig. 5). This is due to the same reason as 
above: Coma cluster has higher temperature and the data, 
collected from a large field, has bigger statistical significance, 
hence smaller errors. 



4 RESULTS 

In this Section we summarize the restrictions on parameters 
of sterile neutrino, obtained from analysis of the Virgo and 
Coma cluster of galaxies. Based on analysis of the core and 
peripheral regions of these clusters we obtain the following 
results: 

(i) The most robust restriction is based on the require- 
ment that the flux of DM did not exceed the total X-ray 
flux in the given energy bin (chosen to be AE — 200 eV in 
our case). The best restriction comes from the observation 
of the Coma periphery (XMM-Newton observation "Coma 
3" (ID 0124710301) with EPIC EPN camera). It provides 
an exclusion plot, shown in the solid red line on the Fig. 8. 

(ii) For rris > 14 keV, the restriction from the Virgo clus- 
ter becomes stronger (due to the fact, that the temperature 
of the intracluster gas in the Virgo cluster is several times 
lower, than that of the Coma. This results in an improve- 
ment, shown in the green dashed line, continuing previous 
bound above rris — 14 keV on Fig. 8. 



(iii) By analyzing the form of the spectrum, one can re- 
flne the restrictions on parameters of sterile neutrino (as 
discussed in details in Section 3.2). From all the data used 
in this work, the best result was obtained in analysis of the 
Coma center region (see Section 3.2.1 and Appendix A). 
The resulting exclusion region is shown as dashed blue line 
on the Fig. 8. It provides an improvement by a factor 2-4 
over the XRB result of Boyarsky ot al. (200.5) in the region 
2 keV < rUs < 10 keV. One can qualitatively understand 
why Coma center observations gives the best result for sta- 
tistical analysis. As discussed in the Section 3.2.1, the more 
monotonic is the data, the better the method of statistical 
analysis works. The result significantly improves over the 
total flux restriction only in the regions where "there are no 
lines". Due to the fact that Coma temperature is 4 times 
bigger than that of Virgo, the phenomenological model for 
the Coma cluster requires only three additional lines (as dis- 
cussed in details in Appendix A). 



5 DISCUSSION. 

In this work we analyzed the restrictions on the parameters 
of the sterile neutrino, as a dark matter candidate, coming 
from the observations of several clusters of galaxies. Below 
we compare our results with those of other authors and dis- 
cuss what improvements on restriction plots one expects to 
obtain from galaxy clusters data. 

(i) We see that this work brings improvement by a factor 
2-4 (in the region 2 keV < rUs < 10 keV) over the previous 
analysis of Boyarsky ct al. (2005) (see Fig. 8). This should 
be compared with the results of Abazajian ct al. (2001); 
Abazajian (2005), where stronger constraints, derived from 
the restriction of the total flux of the Virgo cluster, were 
claimed. Below, we explain the reason for this difference. 
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Assuming ~ 6 x 10~^, Abazajian ct al. (2001) ar- 
gued that if the sterile neutrino mass would be > 5 keV, 
the "enormous" DM decay line observed in the core of 
Virgo galaxy cluster (or, more precisely from its cD galaxy 
M87) would dominate over the continuum flux. To illustrate 
this, Abazajian ct al. (2001) modeled (via WEBSPEC inter- 
face) the emission of the M87 by MEKAL model with the 
temperature T = 2.5 keV and overall normalization of the 
flux in the 2-10 keV interval equal to 1.5 x 10~^^erg/ (cm^ -s), 
the number taken from Bochringor et al. (2001). This value 
for the emission flux from the center of Virgo is, how- 
ever, two orders of magnitude below the actual flux, as one 
can see from the literature (Fabricant & Gorcnstoin 1983; 
Arnaud t Evrard 1999) or by comparing these numbers 
with the actual experimental data (e.g. observation of M87 
with XMM in 2001 (obs. ID 0114120101). From the latter 
one can measure that the flux in the range 2-10 keV is some 
two orders of magnitude higher (~ 10"^*^ erg/cm^ ■ s). While 
the discussed in the paper of Abazajian ct al. (2001) flux of 
DM is taken from the whole FoV of Chandra (which has an 
area ~ 300 arcmin^), the number 1.5 x lO"^'^ erg/cm'^ • s 
in the Bochringcr ct al. (2001) refers to the emission from 
the active nucleus of the central galaxy of Virgo cluster, 
M87, collected from the central 10 arcsec of the XMM (not 
Chandra) field of view. For the mixing 6) ~ 6 X 10~^ used 
in Abazajian ct al. (2001), the flux of DM is about two or- 
ders of magnitude below the MEKAL flux in the energy bin 
around E = ma/2 = 2.5 keV (see Fig. 9). In the later work 
Abazajian (2005), the same data was used to provide an ex- 
clusion plot in the (nis, sin^ 26) parameter space of the form 
sin2(26l) < 1.0 X 10"^( keV/m,)^ 

(ii) Statistical analysis of the Coma center and Coma 3 
data, gave very close restrictions. However, Coma 3 repre- 
sents only one of the observations of Coma mosaic, which 
cover the annular region between ~ 10' and 40' around the 
center of the Coma cluster. By combining these results, one 
effectively increases observation time by about a factor of 
10. This should lead to the modest improvement of the sta- 
tistical restriction by about \/To. We leave this for the future 
work. 

(iii) The restrictions, based on the total flux, can be im- 
proved if one moves to the periphery of the cluster. The far- 
ther away from the center one moves, the stronger becomes 
"signal-to-background" ratio (especially in the case of the 
Coma cluster, c.f. Eq. (18)). Therefore, one should try to 
move away from the center of the cluster. The maximal dis- 
tance where one can move i restricted by the requirement 
that the DM decay line is still stronger, than the XRB back- 
ground (otherwise one does not gain anything, compared to 
the restriction (4) from XRB). 

(iv) At the energies above 15 keV, other missions (such 
as INTEGRAL) should be used to obtain the restrictions on 
parameters of sterile neutrino. 
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APPENDIX A: SIMPLE MODEL FOR THE 
X-RAY SPECTRUM OF THE CORE OF COMA 
CLUSTER. 

The simple phenomenological flt to the XMM-Newton PN 
and Beppo-SAX PDS data is shown in Fig. Al. The model 
consists of a powerlaw with a photon index F — 1.36 ± 0.04 
modified at high energies by an exponential cutoff with the 
cut-off energy Ecut ~ 8.3 ± 0.7 keV and by the absorption 
with column density Nh = (1.5 ± 0.3) x 10^" cm"^ (which 
is roughly the galactic absorption in the direction of Coma 
cluster). Several bright X-ray lines are visible in addition to 
the cut-off powerlaw component. They are the red-shifted 
6.7 keV and 7.0 keV Fe Ka and K/J lines as weU as 0.5 keV 
and 1.0 keV O and Ne lines. All the components of the model 
are shown together with the unfolded 0.3-50 keV spectrum 
of the central part of the Coma cluster in the right panel of 
Fig. Al. The reduced of the fit is = 1-04. 



APPENDIX B: SIMPLE MODEL FOR THE 
X-RAY SPECTRUM OF THE CORE OF VIRGO 
CLUSTER. 

Modeling the spectrum of the central part of Virgo clus- 
ter turns out to be more complicated because more lines 
are evident in the 1-5 keV part of the spectrum. This can 
be explained by lower, compared to Coma, temeperature of 
Virgo cluster, T ~ 2 — 3 keV. Besides, Virgo cluster is a 
"cD galaxy dominated" cluster, which means that the dif- 
fuse X-ray emission from the cluster is strongly dominated 
by emission from the central giant elliptical galaxy M87. 
The size of the X-ray core of M87 is much smaller than the 
size of X-ray core of Coma cluster, so that in spite of the 
fact that Virgo is about 5 times closer to the Earth than 
Coma, the angular size of its core, ro — 2 arcmin, is 5 times 
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Figure Al. Folded (a) and unfolded (b) XMM-Newton and Beppo-SAX spectrum of the Coma galaxy cluster. The simple phenomeno- 
logical model used to fit the data consists of a cutoff powerlaw modified by the galactic absorption and four Gaussian lines. The model 
parameter values are given in the text. 



less than the one for Coma. This, in turn, implies that the 
temperature of the intracluster (or intragalactic) medium 
varies significantly over the XMM FoV, which extends up to 
6 — 7ro. Finally, Virgo cluster is not a relaxed galaxy clus- 
ter, and assumptions about hydrostatic equilibrium, which 
one has to evoque to compute the DM distribution may not 
hold. 

The above facts make the modelling of the X-ray spec- 
trum collected from the entire XMM FoV difficult. However, 
we achieve a satisfactory phenomenological fit the the Virgo 
cluster spectrum with the model similar to the one used in 
the previous Section for Coma cluster analysis. It consists 
of the powerlaw with the photon index F — 1.28 modified at 
high energies by an exponential cut-off at Eq = 2.22 keV and 
at low energies by photoelectric absorption on the Galactic 
hydrogen column density Nh = 2.5 x 10^°cm~^. Several 
lines have to be added to the fit, to make it acceptable. In 
particular, Fe Ka and K/3 lines of Fe (6.7 and 7.9 keV), line 
of O at 0.64 keV, a set of lines at 1 keV from the cooling 
flow in the center of Virgo cluster, etc. (7 lines altogether). 



APPENDIX C: SIMPLE MODEL FOR THE 
X-RAY SPECTRUM OF THE RING WITH THE 
RADIUS 9' 11' AROUND M87. 

We repeat the procedure already used for the analysis of 
the central parts of Coma and Virgo clusters for the present 
case. Namely, we find a simple phenomenological fit to the 
broad-band spectrum which in this case is given by the cutoff 
powerlaw with F = 1.29 and Eo = 2.24, modified by the 
Galactic absorption. Again, we add several lines (as in the 
case of the previous Appendix) to make fit acceptable. 
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